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Thermodynamic properties of cubic Heisenberg ferromagnets with competing exchange interac- 
tions are considered near the frustration point where the coefficient p in the spin-wave spectrum 
-Ek ~ pk 2 vanishes. Within the Dyson-Maleev formalism it is found that at low temperatures 
thermal fluctuations stabilize ferromagnetism by increasing the value of p. For not too strong frus- 
tration this leads to an unusual "concave" shape of the temperature dependence of magnetization, 
which is in agreement with experimental data on the europium chalcogenides. The phase diagram 
is constructed by means of Monte Carlo simulation, and suppression of magnetization and Curie 
temperature is found in comparison with the results of the spin- wave theory. This effect is explained 
by the the presence of non-analytical corrections to the spin-wave spectrum which are represented 
in the lowest order by the term ~ (T/S) 2 k 2 logfc. 
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During last thirty years, the physics of quantum mag- 
nets has been developed under strong influence of the 
concept of quantum disordered ground state, which gives 
an alternative to the standard magnetic order picture. 
Main results in this direction were obtained for antifer- 
romagnets in one and two dimensions (see, e.g., Ref. 
1). However, recently the existence of a quantum dis- 
ordered ground state in two-dimensional magnets with 
predominantly ferromagnetic interactions was demon- 
strated, see Refs. and The formation of the cor- 
responding quantum state, called a spin nematic, is as- 
sociated with an instability with respect to the creation 
of bound states of spin waves. Since in the ferromag- 
netic ground state quantum fluctuations are completely 
absent, the physics of the "frustrated ferromagnets" dif- 
fers substantially from the physics of antiferromagnets 
with competing exchange interactions. 

Magnetic frustrations also significantly influence the 
thermodynamic properties of three-dimensional systems, 
in particular the temperature dependence of long- and 
short-range order parameters (see e.g. Ref. for the fee 
antiferromagnet). Anomalous temperature dependences 
of the magnetization are indeed observed experimentally 
in ferromagnetic materials, for example in the overdoped 
europium chalcogenides, see Ref. Typically, these 
anomalies are explained by introducing into the Hamilto- 
nian four-spin (biquadratic) interactions^. In this paper 
we show that these anomalies can be naturally obtained 
in the framework of the standard Heisenberg model by 
taking into account the frustration consistently. 

We consider the Heisenberg model to investigate ther- 
mal properties of three-dimensional frustrated ferromag- 
nets. The frustration is caused by the exchange interac- 
tions in higher coordination spheres. Low-temperature 
properties of this model can be determined not only by 
the conventional magnons, but also by the bound states 
of spin waves 7 ^. However, as it was shown in Ref. 0, in 



three dimensions frustration does not lead to the forma- 
tion of the bound state of two spin waves at long wave- 
lengths. Therefore, the contribution of the bound states 
is absent at low temperatures. 

The spin- wave spectrum Eo(k) = S [ J(0) — J(k)] 
(known exactly for a ferromagnet at zero temperature) is 
quadratic in the wave vector at k — »■ 0, EqQs.) s» p a pk a kp, 
where p a p — S^2 r J r r a rp/Q is the spin- wave stiffness 
tensor. Here J(k) is the Fourier transform of the ex- 
change interaction J r , r a are the coordinates of the lat- 
tice sites. Provided that owing to the frustration two 
or three eigenvalues of p a p vanish, the correction to the 
magnetization in the spin-wave theory, 



Mswt = S — 



d 3 k 



1 



(2tt)3 exp [£ (k)/T] -I 



(1) 



diverges at any finite temperature T. The exact form of 
the divergence is determined by the higher order terms 
of the expansion of Eo(k) in k. Obviously, when all three 
eigenvalues of p a p vanish simultaneously, the divergence 
is stronger. This case is easily realized for cubic crystals 
where the tensor p a p = po 8 a p is always diagonal. Below 
we assume the presence of cubic symmetry. 

Corrections to the spin- wave theory can be calculated 
systematically using the Dyson-Maleev representation 
for spin operators^. The corresponding effective boson 
Hamiltonian is H c jf = Hq + H mt where 



Hv = Y^E Q (y)a{a v 



(2) 



is the Hamiltonian of non-interacting magnons, 

Hint = j ^2 !/>(!, 2; 3,4) 44a 3 a 4 <5i +2 ,3+4> ( 3 ) 



1,2,3,4 



is the interaction, (p (1, 2; 3, 4) = J (3) + J(4) - J(l - 3) - 
J(l —4). The spin- wave spectrum at finite temperatures, 
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E = -Eo+Re E, is determined by the self-energy of bosons 
E(k, iui n ). Fig. 1 shows low-order Feynman diagrams for 
E. 




a b 



FIG. 1. Diagrams for boson self-energy in the first and second 
orders in 

Let us consider the contribution of the diagram (a) to 
the spin- wave stiffness at finite temperature, p(T). Direct 
evaluation shows that at low T 

p(T^ Po )=p Q -A(T/po) 5/2 . (4) 

Under normal conditions, i.e. in the absence of the frus- 
tration, the coefficient A > 0, so that with increasing 
temperature the spin-wave stiffness decreases and the 
magnetic order is destroyed. However, if the frustration 
is present, the situation may chang o 9 ' 10 . Note that the 
coefficient A can be expressed in terms of the coefficients 
of the expansion of E$ (k) in k. For cubic crystals there 
are only two fourth-order invariants, 

Ea(k) « Po k 2 + *ik 4 + x 2 (fc^ + k* + fc*). (5) 

Then A = -C(5/2)(5xi + 3x 2 )/87r 3 / 2 S. We assume be- 
low that even for an arbitrarily small po the ferromag- 
netic ground state remains stable. The requirement of 
the positivity of Eq. (J5J) for po — > then gives inequal- 
ities x\ + k% > 0, 3xi + ?*2 > 0. One can see that in 
the considering case these inequalities lead to A < 0. For 
po = we find a behavior, analogous to Eq. ((4]), but with 
a different exponent, 

P (T « x) = -bA (T/x) 5/4 (6) 

where x = max{xi,X2}, b > is the function of xi : 2- 
Rearranging the diagram series in the occupation num- 
bers of bosonsii, one can show that the dependences ((4]) 
and (JSJ) do not change their form in higher orders in 1/S. 
Thus, regardless of the details of the exchange interac- 
tions, in the frustration regime the ferromagnetic order 
at low temperatures is stabilized by thermal fluctuations 
(similarly to the "order from disorder" effect, see, e.g., 
Ref. H2). 

The nonlinear spin-wave theory (which takes into ac- 
count only the diagram Fig. la) is valid in the limit 
T — > 0. We have furthermore performed the calculations 
within the self-consistent spin- wave theory (SSWT), pre- 
viously used for two- and three-dimensional systems^. 
SSWT includes all kinds of diagrams obtained from the 
diagram (a) by its recursive reinserting into internal 
lines. To improve the behavior of SSWT results at high 
temperatures, the contributions from pseudofermions are 



included^, the chemical potential of bosons being intro- 
duced to fulfil the equation M = in the paramagnetic 
phase 15 . Fig. 2 shows an example of dependences p(T) 
calculated in SSWT for 5=1/2 Heisenbcrg model on the 
fee cubic lattice. It can be seen that at low temperatures 
near the frustration point po = the dependences p(T) 
are consistent with Eqs. (0J and ©. 
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FIG. 2. Temperature dependence of the spin-wave stiff- 
ness calculated in SSWT for S — 1/2 Heisenberg model 
on the fee cubic lattice with the interaction between the 
first, second, and third neighbors. J2 = — 0.3Ji, and J3 = 
(po /S — Ji — J2 ) /6 was chosen for a number of values of po 
shown in the figure. The filled circles correspond to the Curie 
temperatures. 

Note that the spectrum E(k) in SSWT equation for 
the magnetization, Eq. ([1}, is renormalized by temper- 
ature corrections. As a result, the magnetization does 
not diverge at po = 0. However, the growth of the spin- 
wave stiffness with increasing temperature, see Eq. ([6]), 
results in this case in unusual behavior of the magneti- 
zation, M ~ S — uT 3 / 8 . Fig. 3 shows the temperature 
dependence of the magnetization for several values of po. 
It can be seen that for small non-negative po the curves 
M(T) at low temperatures acquire a concave form. Such 
anomalous dependences were observed experimentally in 
overdoped EuO^. 

The results for spin-wave stiffness and magnetization 
at negative po (when the ferromagnetic ground state be- 
comes unstable) are also shown in Figs. 2 and 3. In 
this case the temperature corrections lead to the stabi- 
lization of the ferromagnetic ordered phase above some 
"lower" critical temperature. Therefore, the dependence 
of the Curie temperature T c on po becomes two- valued at 
po < 0, see Fig. 4 for a classical system (in the quantum 
case the dependence T c (po) has a similar form). 

When p < and ferromagnetism is unstable, another 
state, which is able to compete with the ferromagnetic 
state, should appear. In principle, this could be a spin 
liquid state. However, the most likely candidate is a spi- 
ral state. Note that for p < the spin-wave spectrum 
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FIG. 3. The dependence of magnetization on temperature in 
SSWT. The model parameters are chosen as described in the 
caption to Fig. 2. 



([5]) has a minimum at k = Q = (±Q, ±<3, ±<9) where 

Q = (-p/2[3xi + ^ 2 ]) 1/2 - This could mean the forma- 
tion of a spiral state with the wave vector Q. This state 
appears to be favourable in particular in the mean field 
approximation. Note that the ferromagnetic and para- 
magnetic states predicted in SSWT without account for 
spirals (which account is difficult) can be thermodynami- 
cally unstable. It is important that the spin-wave theory 
allows to exclude the continuous transition from ferro- 
magnetic to spiral state because this transition would re- 
quire p(T) to vanish at the boundary between two phases. 
The latter is impossible due to the divergence in the mag- 
netization, Eq. (HJ, for p(T) -4 and T > 0. 

For a complete study of the phase diagram and deter- 
mination of the role of non-spin-wave contributions, we 
have performed Monte Carlo simulation of the classical 
Heisenberg model on the finite fee cubic lattices of extent 
L = 10, 14 in each direction with the help of the ALPS 
libraryiS. The resulting phase diagram is shown in Fig. 
4. It contains ferromagnetic (F), spiral (S), and param- 
agnetic (P) phases. The point of frustration po = pf(L), 
which defines the boundary between the ferromagnetic 
and the spiral state at T = 0, is slightly shifted to nega- 
tive po due to the discreteness of the quasimomcntum of 
finite systems. We have found that, for po < Pf and at 
temperatures below the ferromagnetic phase, the para- 
magnetic phase does not arise. Instead, the first-order 
transition to the spiral phase occurs (we have observed 
a two-peak structure of the probability distribution for 
the magnetization). This means that for a given finite 
system the curve T c (po,L) near T = lies in the region 
of thermodynamically unstable states, i.e. below and to 
the left from the line of first-order phase transitions. At 
the same time, in the quantum case the curve T c (p) near 
T = may correspond to the transition to the quantum 
disordered phase. 



FIG. 4. Phase diagram of the classical Heisenberg model on 
the fee cubic lattice with the interaction of the first, second, 
and third neighbors, J2 = — 0.3Ji, J3 = (po — Ji — J2) /6. 
Solid line is the result of Monte Carlo simulation for 10 x 
10 x 10 system. Dotted line is the mean field approximation. 
Dashed line shows the dependence of the Curie temperature 
on po in SSWT. The shaded area is the part of the ferromag- 
netic phase with p(T) < O.Ui according to SSWT. 



For comparison, Fig. 4 also shows the phase bound- 
aries in the mean field approximation. A peculiarity of 
this approximation is the existence of the multicritical 
Lifshitz points which is the endpoint of the boundary 
between ferromagnetic and paramagnetic phases. In the 
Ginzburg-Landau expansion at this point the coefficient 
c oc p(T) at [VM(i)] 2 vanishes. This would lead to a 
modification in the critical behavior* 7 . However, tak- 
ing into account fluctuation corrections to the mean field 
theory in SSWT makes impossible the existence of the 
Lifshitz point at finite temperatures due to the diver- 
gence in Eq. (JTJ for p{T) -> 0. On the curve T c (p ) 
describing the boundary between the ferromagnetic and 
paramagnetic phases in SSWT, the coefficient p(T) is fi- 
nite everywhere and vanishes only in the limit T — > 0. In 
this sense, one can say that fluctuation corrections shift 
the Lifshitz point along the curve T c (po) to the point 
with coordinates T = +0, po = — 0. Infinitely small de- 
viation of T and po from zero is essential, since strictly 
at T — even in the quantum case the system possesses 
saturated ferromagnetism for all po > 0. Thereby we call 
the point with coordinates T = +0, po = — the quasi- 
Lifshitz point. In SSWT, when moving away from the 
quasi-Lifshitz point along the curve T c (po), the size of 
the corresponding "critical" temperature range (where 
the magnetization changes from its maximum value to 
zero) increases as (— po) 4 ^ 5 (see also Fig. 3). There- 
fore, given that the curve T c (po) in the neighborhood of 
T = falls into the region of thermodynamically unsta- 
ble states, fluctuations in the region of stable states are 
expected to be important at sufficient distance from the 
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quasi-Lifshitz point near the curve T c (po). The quasi- 
Lifshitz point affects also the behavior of the magnetiza- 
tion for weakly positive po, see Fig. 3. 




FIG. 5. Temperature dependence of the magnetization for 
the classical Heisenberg model on the fee cubic lattice of fi- 
nite size 14 x 14 x 14. The model parameters are chosen as 
in the caption to Fig. 4. The inset shows temperature depen- 
dences of the short-range order parameter. Solid lines are the 
results of Monte Carlo simulation. Dashed lines are the result 
of SSWT. The filled circles indicate the Curie temperatures 
determined by the maximums of the susceptibility. 

Everywhere, except for sufficiently large po m Ji, the 
Curie temperature in SSWT and mean-field approxima- 
tion is overestimated by 2-3 times in comparison with 
the results of the numerical simulations, see Fig. 4 (cf. 
Ref. M). Fig. 5 shows that for p w J x SSWT gives a 
good description of the temperature dependence of the 
magnetization even at high temperatures, except for a 
narrow critical region near T c . In the presence of frus- 
tration for pf(L) < po <C Ji the ferromagnetic order in 
the ground state is stable, and SSWT still reliable at low 
temperatures, but becomes inapplicable at higher tem- 
peratures. The latter fact is manifested, in particular, in 
a significant overestimation of the Curie temperature. In 
the presence of strong frustration for p < Pf(L) ferro- 
magnetism in Monte Carlo simulation is observed only 
at intermediate temperatures, which is in a qualitative 
agreement with SSWT, but the quantitative discrepan- 
cies appear now at all temperatures. 

The inset in Fig. 5 shows the temperature dependence 
of the short-range order parameter 7 = ^/S.; • Sj (i and 
j are nearest neighbors). In the presence of frustration 
the short range order persists up to T c and above, as 
one would expect. It is remarkable that even in the case 
of strong frustration and high temperatures SSWT de- 
scribes the short-range order much better than the long- 
range one. This means that the basic defect of SSWT 



for strong frustration is associated with insufficient ac- 
count for long- wavelength fluctuations. The existence of 
short-range order can be crucial for the energy gain of 
certain phases, which, in particular, is of interest for the 
description of the j-a transition in ironic. 

Note that when the system is frustrated, i.e. for 
po "C Jx, according to the results of Fig. 2 the spin- 
wave stiffness p(T) is small even at temperatures up to 
the Curie temperature and above. This can explain the 
failure of SSWT in the case of strong frustration. If we 
formally consider p(T) as an independent parameter, a 
perturbation series in 1/5* for various physical quantities 
will contain divergent terms in the limit p(T) — > 0. An 
example is given by the Eq. (Q~]) for the magnetization. 
The singularities appear also in the expansion of the self- 
energy. Consider, e.g., two- loop diagram (Fig. lb). In 
the limit p(T) — >• it gives singular term in the spin-wave 
spectrum 

11T 2 k 

£(v^M « k « A) ~ fc 4 + Q67r2S2xi k 2 log - (7) 

(x2 = 0, A is an ultraviolet cutoff; the condition for 
the applicability of the perturbation theory gives an ad- 
ditional estimate for the momentum, fc 2 /log(A/fc) 3> 
0.01(T/5xi) 2 ). Note that the logarithmic term is neg- 
ative. Similar logarithmic terms are typical for the sys- 
tems with strong fluctuations (see, e.g., Refs. HH and 

El). 

It is important that the singular terms similar to that 
in Eq. ([7]) can result in a significant suppression of the 
magnetic order even at finite but small p(T). Suppression 
of the magnetic order should be particularly strong near 
the lower critical temperature (see the turn of the curve 
T c (po) in Fig. 4) since the value of p{T) is smaller in this 
case (see Fig. 2). 

The results, similar to those discussed in the paper, can 
be obtained for all cubic lattices (simple cubic, bec, and 
fee) with different types of competing interactions. The 
presence of low-temperature regions with strong fluctu- 
ations on the phase diagram could favor the formation 
of a quantum disordered state in three-dimensional frus- 
trated ferromagnets. To achieve advance in this direc- 
tion, taking into account quantum effects in the spiral 
state near the boundary with ferromagnetic phase is of 
considerable interest. The "bold" diagrammatic Monte 
Carlo technique developed quite recently in Ref. can 
provide the appropriate tool for this. 
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